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1. INTRODUCTION 

Recent explosion on the Internet growth has put pressure on optical fiber network to expand its 
channels. To fulfill the requirements on increasing information carrying capacity of optical fiber 
communication systems the dense wavelength division multiplexing (DWDM) technology might be utilized. In 
DWDM system a number of devices have been used which include laser sources with different optical 
wavelength, multiplexer/demultiplexer to combine and spread the signals at the transmitter and receiver side. 
DWDM technology also offers the flexibility to enhance network capacity without adding more fiber optics 
into the existing systems. Therefore, the DWDM can further increase the channel data rates or increasing 
the number of multiplexed channels [1]. 

In principle, DWDM devices that can be used as a set of optical channels, each of which uses 
a different wavelength [2]. Interleavers are becoming an interesting research focus on the DWD systems 
because it can design to provide new effective solutions. An optical interleaver is a device that combines two 
DWDM channels into multiple signals stream in interleaving ways or vice versa [3]. Interleaver light signals 
work similar to comb filters, it can split one signal into two or more signals with greater intervals between 
them [4]. A method is often used to design filters in optical fiber is the Interleaver method based on Mach-Zehnder 
interferometer (MZIJ) structure [5]. 
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The MZI has been developed to act as a switch (in its balanced form) and as a multiplexer or 
demultiplexer (in its unbalanced form). The MZI consists of two directional couplers or Y-junctions, 
connected by linking fiber optics. Introducing a phase delay to one linking guide via the thermo-optic effect 
enables the MZI to be used as a switching device. Based on the number of couplers, the configuration of MZI 
is divided into several types, namely single MZI consists of two couplers and cascaded MZI use more than 
two couplers [5-7]. The characteristics and the design of MZI have been the focus of discussion in 
developing sensitive filters and sensors [8, 9]. 

In this paper, cascaded MZI structure is used as optical interleaver with one input, namely 
on port | [9, 10]. In general, interleaver devices can operate well in a certain temperature range. However, 
outside the predetermined temperature range, the signal transmission characteristics on the interleaver device 
will change [11, 12]. Based on this problem, we simulated the effect of heating two-stage and three-stage 
successive MZI with the range of temperature from 30°C to 430°C. 


2. RESEARCH METHOD 
2.1. Temperatures effects on silica fibers 

The thermo-optics (TO) effect is refractive index changes inside the fiber optics as the temperature 
rises [12, 13]. Usually, the TO phase shifters are comprises of a thin film heater placed on the fiber optics. 
When the heater is on the material refractive index will change, for example heating a 10 mm long guide at 
temperature of 15.23°C, will create a z radians phase shift at 1.523 um wavelength [8]. The phase shift 
values and time to heat the materials depend on the cladding thickness, thermal characteristics, and substrate 
material used [9]. 

At this time, there has been a new research to determine the new relation in equation 
and thermal-optical coefficients in silica fibers at temperatures ranges of 0-1200°C. The relation also includes 
second order coefficients of thermo-optics as follow [14]. 


n(T)=ng+a,T+B T° (1) 


where n(T) is the effective refractive index as a function of temperature. a, and f, are first-order 
and second-order thermo-optic coefficients respectively. In this calculation we adopt the values 


as a,=1.090x 10° and p =1.611* 10° [14]. Figure 1 shows the effect of temperature on the effective 
refractive index in optical fiber. The changes in refractive index are a linear to the increase in temperature. 
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Figure 1. The changes of effective refractive index as temperatures increase 


2.2. Fiber optics Mach-Zehnder interferometers 

The structure of MZI is generally divided into two parts, namely symmetrical MZI and asymmetric 
MZI [12, 15, 16]. Symmetric MZI has the same length of the upper and lower arm, whereas in MZI 
the asymmetrical arms have different lengths. In this paper, we just focus on asymmetric MZI. MZI consists 
of two or more 2x2 directional couplers which are connected by different lengths in each arm [17-18]. 
The difference of length made a light wave in one of the optical fiber paths delayed, thus creating a phase 
shift on the two light beams at the output [19-21]. In this section, MZI consists of two couplers that 
are connected together on their both ends using optical fiber. This structure is referred to as a single MZI [6, 7]. 
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Figure 2 shows the structure of MZI with two couplers, the input signal is launched into port 1 or 
port 2 and the output can be detected at point 3 and 4. In this case L is the length of MZI arms and AL 
is the different length of both arms. The C; is coupling ratio of the first coupler and C; is the coupling ratio of 
the second coupler. The input light will be divided into two output at first coupler and travels in a different 
arm length that create a phase sift between the two modes. Further down the two modes will be recombine 
by coupler 2 and split at the output arms [9, 16]. 







L+AL 
Port 1 Port 3 
Input Output 
Port 2 L Port 4 


Figure 2. Two couplers MZI diagrams (single MZD) 


The output or transmission characteristics of MZI can be explained using the transfer matrix 
method [16]. It is divided into two main matrices namely couplers and phase shift matrices. The couplers 
transfer matrix is obtained from [6]: 


cos -jsin @ 
i | (2) 


me=| 
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where 0=KL,, K is coupling coefficient in the coupler and Le is coupling length in the coupler. 

This parameter is used to determine the value of the coupling ratio used in the MZI structure [21-23]. 

Meanwhile, a phase shift matrix can be represented as follows [16, 18]: 
_fe* 0 

“=| 0 et (3) 


where [11]: 
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where ø, is phase shifting on single MZI and £ is the mode propagation constant in the delay section that 
determined by: 


p= ney (5) 


Here nep is effective refractive index of optical fiber. Transfer matrixes for single MZI structure are 
consisting of two transfer matrixes. First is transfer matrixes with coupling coefficient on two couplers and 
the second is transfer matrix of phase shift (Mọ) as shown in Figure 2. So that it can be written as follows: 


m=|® is fen 0 | |2 A (6) 
JS] C] 0 ei) |-js2 c 


where, s;=sin 0; =./C,;, c;=cos 6; = 11-C,, and i=1,2,3...n (number of couplers) and C; is a value 


of coupling ratio in port 3. The equation above is the transfer matrix equation for single MZI structure or 
MZI structure with two couplers [20]. To get an output signal the following calculation might be used as: 
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The powers outputs for port 3 and port 4 can be obtained by calculate (8) yielding: 
Poutt=lEoutt |? and Powt2=\Eout2|? (9) 


Optical Interleaver is a device that can be scaled, so its spectral properties are periodic [9, 12]. 
The periodicity of the interleaver is referred to as the free spectral range (FSR) 1.e. the optical frequency 
distance or length of two waves of optical intensity that are reflected or transmitted sequentially from 
an interferometer. Separation of a standard ITU 50 GHz channel requires 29=1550 nm. The relationship 
between differences in arm length and FSR 1s [12]: 


32 
FSR = — 
Neff ÁL 





(10) 


where /, is the wavelength center, nepis the effective refractive index at that particular wavelength, and AL is 
the arm different length between two arms of MZI [22, 23]. Hence, (10) might be used to get the length 
difference needed to obtain a specific FSR at a specific wavelength. 


3. RESULTS AND DISCUSSION 

In this calculation we will analyze several MZI structures that consist of one stage MZI, two-stage 
cascaded MZI structures, and three-stage cascaded MZI. The output power of each MZI structure are 
calculated similar to those of single stage MZI except that for each additional cascaded MZI we have to insert 
a new coupler and phase shift matrix. Furthermore, we also analyze the output characteristics of those MZI 
structures as a function of temperature changes. 


3.1. Temperature effects on two-stage cascaded MZI 


Figure 3 shows the structure of the two-stages cascaded MZI [9] which consists of three directional 
couplers and second phase shifters where the second delay line is 24Z. 


Port 1 Port 3 






Input Output 


Port 2 Port 4 


Figure 3. Two-stage cascaded MZI structure with heating element on arms 


Therefore, the output field from (4) can be revises to become: 


_ B2AL 


$= (11) 


Using the similar method as (6) the output power of the two-stage Cascaded MZI is derive as: 


wje He sje +e she ay (2) 
-JS] C] 0 ot -JS2 C? 0 eJ? -JS3 C3 


To simulate the power output of the two-stage cascade of the MZI we use the following parameters 
as shown in Table 1. The coupler coupling ratio is determined based on coupling coefficient and coupling 
length in 2x2 optical fiber couplers [16]. To get the spectral response on port 3 and port 4 of two-stage 
cascade MZI, we have use (9, 12). Figure 4 demonstrates the spectral response of two-stages cascaded MZI 
with the C-Band wavelength spectrum. Spectral responses on the two ports have similar characteristics, 
where the signal is isolated in the same wavelength range at each particular wavelength. This isolated signal 
range is called the isolation channel. Separation of channels in each isolation channel is called a crosstalk 
channel. Directivity or crosstalk is a measure of power reflected from an MZI device. This explains that the 
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coupler is very good for transmitting light to the output port and calculating the amount of light reflected in 
the device [18]. When the crosstalk value is very high it will interfere with propagation. The characteristics 
of Figure 4 are summarized in Table 2. 

One way of changing the silica fiber refractive index is by using the thermo-optical effect. 
These changes make the loss of transmission low, but still retain its basic features [4]. The temperature effect 
can also be used as a sensor that has high sensitivity, fast response, and immunity to electromagnetic 
interference [15]. The effective thermo-optical coefficient depends on the fiber material, temperature changes 
in each layer, and fiber diameter [14]. The transmission power characteristics of two-stages cascaded MZI as 
a function of temperature are simulated with a temperature range between 30°C (room temperature) to 430°C 
with A7=100°C. The output signals are calculated from port 3. Figure 5 shows a shift in the wavelength 
resulting from changes in temperature in port 3. The simulation results in the C-Band wavelength range 
and at the center wavelength of 1550 nm. Heating on the Cascaded MZI Two-stage arm only affects 
the wavelength shift as summarized in Table 2. 


Table 1. The simulation parameters 


Parameters Symbol Unit Value 
Effective refractive index in fiber Neff 1.4645 
C-Band wavelength range A Nm 1530-1565 
Radius of fiber core a um 4.1 
Coupling ratio of 2x2 optical fiber couplers (port 3/port 4) G: % C,=50/50 
C2=75/25 
C3=90/10 
Cz=99/1 
Arm’s length different AL um 178.3 


Table 2 Characteristic of two-stage cascaded MZI 


Result Value 
— Isolation power (port 3 and port 4) -32.60 dB 
a Crosstalk power (port 3 and port 4) -24.50 dB 
i — Pout1 Width of isolation channel 8.20 nm 
$ --- Pout2 Width of crosstalk channel 1.00 nm 
a FSR 9.20 nm 
Center wavelength of FSR 1555 nm 






Space between output at port 3 and port4 4.60 nm 
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Figure 4. Spectral response in port 3 and port 4 of 
two-stage cascaded MZI 
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Figure 5. Spectral response of two-stage cascaded MZI (Port 3) with heating on arms 


3.2. Temperature effects on three stages cascaded MZI 
Figure 6 shows the characteristics of the three-stages cascaded MZI which consists of four couplers 
and three phase shifters where the third delay line is set to be 24L. Similar to the previous analysis, (4) can be 
revised to become: 
aho 


=e (13) 
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Figure 6. Three-stages cascaded MZI structure with heating element on arms 


Using the same method as (6) the output power of the three stages cascaded MZI are derived as follow: 


ele HEE SI HEE Se, E se 24 0 
-jS' C] 0 elt -JS2 C2 0 el, -JS3 C3 0 eT? -JS4 C4 


The spectral response on port 3 and port 4 of two-stage cascade MZI with one input at port 1 
is shown in Figure 7. The results from Figure 7 are summarized in Table 3. The spectral response 
characteristics demonstrate on Table 3 is a bit different from Table 2. The spectral characteristics of 
the three-stage cascade MZI have the isolation power and the crosstalk power lower than that of two-stage 
cascaded MZI, but slightly have a loss of -0.2 dB. Besides that, the width of the crosstalk channel in 
the schematic design looks wider. Output on both ports has a difference in the isolation power values, 
additionally the isolation power value at port 4 is higher than that at port 3. In the same way as in previous 
section and using parameters in Table 1, the transmission power characteristics of three-stages cascaded MZI 
as a function of wavelength are simulated using a temperature range 30°C to 430°C with a range of 
AT=100°C is demonstrated in Figure 8. The simulation based on the wavelength range of 1530 nm to 1565 
nm (C-Band) with the wavelength center at 1550 nm. Value of power isolation, crosstalk power, the width of 
crosstalk and isolation channel, and loss are similar as those in Table 3. Heating on the cascaded MZI with 
three-stage arm can only affects the wavelength shift. 


Table 3. Characteristic of three-stage cascaded MZI 












ee Result Value 
— i wi Isolation power in port 3 and port 4 -37.5 dB and 40.0 dB 
S i-a r i Crosstalk power in port 3 and port 4 -16.0 dB and -16.0 dB 
S ui wy! — Pout 
g u ut --- Pout2 es ae ee 
Š " E Width of isolation channel 7.7 nm 
a. ! i Width of crosstalk channel 1.5 nm 
! i FSR 9.20 nm 
-40 i | Center wavelength of FSR 1550 nm 
1530 1535 1540 1545 1550 1555 1560 1565 Space between output at port 3 and port4 4.60 nm 
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Figure 7. Spectral response in port 3 and port 4 of 
three-stage cascaded MZI 
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Figure 8. Spectral response of two-stage cascaded MZ] at different temperatures 
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3.3. Wavelength shift caused by temperatures 

We have simulated the temperature effect on two-stage and three-stage cascaded MZI 
in the previous discussion. In this section we focus on the temperature effect of two stages and three-stage 
cascaded MZI. Based on simulation result in Figures 5 and 7 we obtained that the heating effect on MZI arms 
that affecting the wavelength shift [23, 24]. It can be seen that the wavelength shift is increased linearly 
with the raise of temperatures. Figure 9 shows that the two stages and three-stage cascaded MZI has similar 
characteristics with the same wavelength shift, therefore we try to find the formulae to get the relation 
of wavelength shift and temperature changes [25]. From the data in Figure 9, we can find the gradient of 
graphics as: 


_ 4836 
~ 430-30 





=0.01209 (15) 


—e— Two-stages CMZI - ©- Three-stage CMZI 


T = 430°C 


4.836 
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Figure 9. Wavelength shift as a function of temperatures 


We can easily obtain the equation between the wavelength shifts due to temperature as [23]: 

Ad=mT + C (16) 
Here m is a gradient, C is a constant. If we substitute (15) into (16), we get: 

(Ar- 29 )=0.01209( Tr - Tọ) (17) 


where Tọ and Tr are respectively room temperature and final temperature, 49 is the wavelength 
at temperatures Tọ and Àr is the wavelength at the final temperature (Tr). As a result we can easily derive 
the formula for wavelength shift due temperatures changes as: 


Ap=Ap + 0.01209( Tr - Ty ) (18) 


Generally, this equation can be used to predict the wavelength shift against temperature on cascaded 
MZI. By using the wavelength shift equation, we can immediately choose the right temperature to be used 
to shift a certain length of wavelengths in cascaded MZI. Further application of this equation it can be used 
to calibrate the wavelength shift so that it may be used as a sensor device based on optical fiber optics. 


4. CONCLUSION 

It can be concluded that the heating effect in two-stages and in three-stages cascaded MZI cause 
the wavelength shift. The shape of the cascaded MZI structure does not affect the value of wavelength shift, 
because the values will be similar on each of structure. Therefore, we can predict wavelength shift by simply 
using the linear variation of wavelength to the temperature changes as shown in. This equation may be used 
as a basis of designing the cascaded MZI based thermo-optic effects so that the characteristics of each 
structure can be predicted in advance. 
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